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Abstract. A noisy damping parameter in the equation of motion of a nonlinear oscillator renders the fixed 
point of the system unstable when the amplitude of the noise is sufficiently large. However, the stability 
diagram of the system can not be predicted from the analysis of the moments of the linearized equation. 
In the case of a white noise, an exact formula for the Lyapunov exponent of the system is derived. We then 
calculate the critical damping for which the nonlinear system becomes unstable. We also characterize the 
intermittent structure of the bifurcated state above threshold and address the effect of temporal correlations 
of the noise by considering an Ornstein-Uhlenbeck noise. 

PACS. 02.50.-r Probability theory, stochastic processes and statistics - 05.40.-a Fluctuation phenomena, 
random process, noise and Brownian motion - 05.10.Gg Stochastic analysis methods (Fokker-Planck, 
Langevin, etc.) 

1 Introduction linearized equation can be misleading : higher moments 

seem to be always unstable although it is known that the 

It is a well-known fact that a multiplicative noise acting -i- i i r ^ i ^ • i-i, i 

^ ° critical value oi the control parameter is the same as that 

on a dynamical system can generate unexpected phenom- deterministic system ^. This apparent contradic- 

ena such as stabilization 1,2 , stochastic transitions and ,.i r-i ^--i-ii ii- 

tion IS due to the existence oi long tails m the stationary 

patterns 3. or stochastic resonance 4 . As far as the sta- i i -r; j- , -i ,• r j-i r • j i- i.- \. 
^ ' ■ ' - - probability distribution oi the linearized equation which 

bihty of a random dynamical system is concerned, the ex- suppressed when the nonlinearities are taken into ac- 

actly solvable example of a nonlinear first-order Langevin ^ rpi , r i ■ j c 

^ count. Ihe same type ot conclusion can be drawn tor a 

equation shows that the behaviour of the moments of the 
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second-order system, namely a nonlinear oscillator with criterion valid for general nonlinear stochastic dynamical 

fluctuating frequency. The energetic stability analysis of systems |12II13| . 

the linearized system has been performed long ago |S1[7] The outline of this work is as follows. In section [3 we 

but it leads to an erroneous phase diagram. In fact, be- define the model and summarize the results obtained from 

cause of the nonlincarities, the noise can even enhance the the stability analysis of the moments. The main part of the 

stable phase m21- Besides, a non-perturbative reentrant paper (section |3Jl is devoted to the calculation of the Lya- 

transition can occur : a noise of small amplitude stabi- punov exponent of the system from which the marginal 

lizes the system, but a strong noise destabilizes it. Similar stability curve is drawn. We then analyze the nature of 

features hold in spatially extended systems such as the the bifurcation and show that the bifurcated state is in- 

Ginzburg-Landau and Swift-Hohenberg equations driven termittent (section 0J|. In the last section, we consider an 

by a noisy control parameter OIS] . Here again, the ob- Ornstein-Uhlenbeck noise to investigate the effect of tem- 

served threshold shift can not be tackled by the analysis poral correlations on the instability threshold, 
of the moments' behaviour alone. 

2 Definition of the model and basic results 

In this paper, we study the nonlinear oscillator with We consider the following stochastic dynamical system: 

random damping. This dynamical system, in which the , 9 ^ , 

^ ^ ^ J ' i + + - i{t)]x + ax = -bx^ , (1) 

noise acts on the velocity variable rather than on the po- 

where ^(i) is a white noise with intensity D and the pa- 

sition, provides another example of a bifurcation induced 

rameters a and h are taken to be positive in order to have 

by a multiplicative noise. This model was introduced in 

a stabilizing effect. The mean value of the damping r is 

the study of the generation of water waves by wind [HI 

also positive whereas the coefficient a of the linear term 

where the turbulent fluctuations in the air flow is mod- 
can be either positive or negative. In fact, as far as the sta- 

eled by a noise. This model is also relevant in the study 

bility of the fixed point is concerned, the a < case can 

of dynamical systems with an advective term where the 

be reduced to the a > case as follows. For a < and in 

corresponding velocity fluctuates. This problem has been 

^^^^ the absence of noise, the oscillation around the x — Q po- 

investigated in |1UII11| via the stability of the moments of 

sition is unstable and the new equilibrium position in the 

the linearized equation. Here, we demonstrate that, as in 

saturated regime is xq = ± \J —ajh. Defining y ^ x — xo, 

the case of the random frequency oscillator, the moments 

we observe that y satisfies: 

stability analysis has no relevance and can not be used to 

derive the instabihty threshold. We shaU obtain the ex- V + {r - + 2axoy + ay^)y - 2ay (2) 

act phase diagram of the system thanks to an instability = x£,{t) — b{3xoy^ — y^) . 
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The linear part of this equation is the same as that of the velocity (whereas the positive fluctuations decrease 

equation once the substitutions r' = r — aa/b and the velocity and have a lesser impact). For the second or- 

a' = —2a are made. The coefficient a' is now positive, der moment, one must consider simultaneously (x^), (i^) 

We thus take a > 0. Before proceeding, we write time and (xi), and study the linear system that couples them. 

and position in dimensionless units by multiplying t and These three quantities become zero in the long time limit 

X by the factors l/\/a and y^ajr, respectively. Defining provided 7 > Z\. fn this case, the system is said to be 'en- 

the following parameters: ergetically' stable. The instability threshold thus depends 

J, jj jjj. on the moment which is considered, due to the interplay 
7 = , A — —= and k = — , (3) 

between noise and nonlinearities. Therefore, in contrast 

equation (|TJ becomes ^j^^ deterministic case, a naive linear stability analysis 

/T , 2\ • , • 7 3 / A \ fails to lead to conclusive results. A more refined criterion 
X + 7(1 + X )x + x = xE^yt) — kx , (4) 

is needed to determine stability. 

where the autocorrelation of the white noise is given by 
{S,it)^{t')) = A6{t — t') (as usual, the brackets represent 

an average on the realizations of the noise). The aim of 3 The Lyapunov exponent 

this work is to study the stability of the fixed point a; = 

_ r 1 1 r 1 1 rr Thc corrcct criterion that determines the stability of the 
of equation IgJ as a function of the values of the dinerent 

fixed point x = is based upon the Lyapunov exponent 

parameters 7, A and k. 

. , , r 1 1 r- 1 defined as 
We now summarize the results tor the stability of the 

moments of the amplitude a;, obtained by Jineariziueequa- , l;o 1 rr 2i\ i'^t(5x) 

yl= -((9fln[fe ]) = (— — ), (6) 

tion ^ around the fixed point x = 0. It can be shown that 

the first moment satisfies : ^^^^^ satisfies equation Q linearized around x 0. It 

2 /\ J been shown |12II13II14| that the sign of the Lyapunov 

-^{x) + ij--)-{x) + {x)^0. (5) 

z ai exponent, calculated with the linear part of the equation, 

The first moment remains bounded as i — s- 00 provided monitors the instability of the nonlinear oscillating sys- 

7 > A/2. The effect of the noise is therefore to enhance tem. When A is negative, the trajectories of the nonlinear 

the unstable phase. An intuitive argument of van Kam- system Q almost surely decay to zero and in the long 

pen |10| explains why this should be the case : positive time limit, the oscillator becomes localized in its rest po- 

and negative fluctuations are equiprobable, but because sition a; = 0. On the contrary, if /I > 0, the flxed point 

the noise is multiplied by the velocity, the negative fluctu- x = is unstable and the stationary probability density 

ations have a stronger effect because they tend to increase of the oscillator is extended. 
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We now calculate the Lyapunov exponent as a function at the origin (the cases z < and z > must be ana- 
of the parameters 7 and A. Defining the variable z as lyzed separately and two different reference points must 



z — dt{5x) / 5x, the equation satisfied by z is given by 

i = -1 - 72 - + . (7) 

This is a Langevin equation involving the variable z alone 
and coupled to a multiphcative noise (to be interpreted 
in the Stratonovich sense). By definition we observe 
that 

A{^,A)^{z). (8) 

Using standard techniques |15II16| . we obtain the Fokker- 
Planck equation for the probability density P{z^ t) at time 
t. The stationary probability density Ps{z) satisfies: 



A 



d,[zPs{2 



1+JZ + Z^ 



Ps{z) = 



(9) 



where J is the constant current of probabihty. Here, the 
current J does not vanish; intuitively, this is due to the 
fact that z plays the role of an angular variable in phase 
space and must be interpreted as a compact variable (see 
[T7j for more details) . Equation @ can be solved by vari- 
ation of constants method and we obtain 



exp 



exp 



A 

2<P{z) 



<P{z) - m) 



dt 



(10) 



where <P{t) = t — 1/t; the constant A and the reference 
point c are not determined at this stage. We must impose 
A — so that Ps{z) is integrable at infinity (this implies 
J 7^ 0). Furthermore, the function exp[f?(t)] is exponen- 
tially divergent when t — s- 0^; this remark fixes the choice 
of the value of c in order to ensure that Ps^z) is integrable 



be chosen). Finally, we obtain. 



J, j:^H{z,t)dt if z<0 



^ H{z,t)dt if z > 



(11) 



where H{z, t) represents the function under the integral 
sign in H10() . Using equation we find that the proba- 
bility density Hill) satisfies 



Ps{z) ^ — for z 



Ps{z) ^ J 



for 



±00 . 
■ 0. 



(12) 
(13) 



This quadratic decay insures that the probability density 
is indeed normalizable at infinity. The normalization con- 
stant N can calculated exactly [ISj and we obtain 



(14) 



■A' 

where J and Y are Bessel functions of the first kind and 
second kind, respectively. This formula was derived in |19| 
for the study of diffusion in a random medium. 

According to equation (jS)), the Lyapunov exponent is 
equal to the mean value of z calculated with respect to 
Ps(^); this quantity will be evaluated in the sense of prin- 
cipal parts 



A = hm 



+M 



zPs{z)dz 



(15) 



M 



(Thanks to equation H12|) . the logarithmic divergencies at 
—CXI and +cxd cancel with each other by parity). An exact 
closed formula for the Lyapunov exponent can then be 
derived and is given by 



A 



NJo 



+ 00 



{l-^)x\dx, (16) 
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where Ki is a modified Bessel function of the second kind 
and the normalization constant N is given in equation \1A\ . 
From this expression, the instability threshold is readily 
found. Indeed, because the Bessel function is always pos- 
itive, the sign of the Lyapunov is the same as that of the 
argument of the hyperbolic sine. We conclude that the 
fixed point x = is stable when 

A 



7 > 



(17) 



We find, as usual, that the stable phase is wider than what 
moments stability predicts. Reverting to the initial vari- 
ables, we observe that the instability threshold is given by 
r — D/4: which does not depend on the 'linear frequency' 
a of the oscillator. In the (r, D) plane, the critical curve in 
the stability diagram of the x = fixed point is simply a 
straight line. We emphasize that this exact result is much 
simpler than that obtained for the Dufhng oscillator with 
random frequency |13j . 

In figure n we plot the Lyapunov exponent of the non- 
linear equation Q for 7 = 1 and 7 = 4. We observe that, 
for yl < 0, the numerical findings and the analytical pre- 
diction are in perfect agreement. For higher values of the 
noise, the nonlinear term can not be neglected : whereas 
the linear Lyapunov exponent, given by equation H16() con- 
tinues to grow, the Lyapunov exponent of the nonlinear 
equation saturates to a value close to 0. This indicates 
that the nonlinear system has gone through a bifurcation : 
the fixed point a; = has become unstable and in the long 
time limit the system reaches an extended stationary state 
(in our case a noisy limit cycle). A remarkable feature of 
the curves displayed in figure ^ is that noise can have a 




Fig. 1. Typical evolution of the Lyapunov exponent as a func- 
tion of the noise amplitude A. The circles/crosses are obtained 
by numerical simulations of the nonlinear equation in the sta- 
tionary regime for 7 = 4 and 7 = 1. The solid lines represent 
the theoretical values llfcil for the linearized equation. 

stabilizing effect at small amplitudes : we notice that the 
curve for 7 = 4 decreases for small values of A and then in- 
creases. A closer inspection shows that this non-moirotonic 
behaviour occurs when 7 > 2. This feature seems to con- 
tradict the argument, presented in section |21 stating that 
a noise in the damping term always has a distabilizing 
effect 10 . This contradiction stems from the asymmetry 
of growth rate of the deterministic system around 7 = 2. 
Indeed, for 7 < 2 the growth rate is a linear function of 7 
((7 = —7/2), whereas it has a very different expression for 



7 > 2 (cr = + — 4/2). The positive and negative 
fluctuations of the growth rate are therefore intrinsically 
asymmetric, in contrast to the implicit assumption in van 
Kampen's argument. 



6 
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Fig. 2. Numerical simulations of the system I^J with 7 = 1, 
k = 1 and A — 3 (top left panel) corresponding to an absorbing 
state and for A — 4.5 (top right panel), corresponding to an 
intermittent state. The lower panels compares the probability 
distribution for z obtained by simulations (dots) with the exact 
result llTTll (full lines). 

4 Intermittency above the threshold 

As shown in figure [3 the temporal series of the sig- 
nal x{t) exhibits "on-ofF intermittency" above the instabil- 
ity threshold (i.e., for a positive Lyapunov exponent). In 
other words, the amplitude x is vanishingiy small for the 
most of the time but exhibits sudden burst of activity. 
This behaviour has already been observed in chaotic sys- 
tems [21)112 1) and in systems driven by multiplicative noise 
j22ll23ll24ll25ll26| . This intermittency can also be identified 
in the probability distribution density. In figure|31 we show 
the probability density of the energy of the oscillator de- 
fined as r = ^/x^ + i^. Near the origin, we observe that 
the probability distribution of r is a power law which is 




r 



Fig. 3. Probability density (in log-log scale) of the variable 
r = \/x^ + for 'y — 1, A — 4.2 and k — 1. The dotted line 
corresponds to the approximation l|18|l . 

generic of an intermittent signal. This power-law distri- 
bution for the energy of the oscillator can be derived as 
follows. From equation Q, we deduce that r satisfies 

r = -{-f + £_{t))r sin^ e - {jam^0cos^ 9 + kcos^ 9 sin e)r^ , 

where 9 — arctan(a;/x) = arctan(z). This equation, to- 
gether with 0, forms a system of stochastic equations. 
A Fokker-Planck equation for the joint density P{r, 9) or 
P(r, z) can be written but the resulting partial differential 
equation seems unlikely to be exactly solvable. 

We therefore make an approximation along the lines 
of |13| and assume that, for small values of r, the proba- 
bility density is separable and can be written as P{r, 9) — 
Ps{r)Ps{9). An average over the angular variable yields an 
independent equation for Ps{r) that can be exactly solved: 
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where the coefficients A, B and C can all be expressed as 
mean values over the variable z (or, equivalently, over the 
angular variable 9) : 

4(1 + Z2)2^ (1 + Z2)2 -^(1 + ^2)2' 

Multiplying equation © by z^/(l + z^) and integrating 
the result over z, the following identity is obtained : 



10 



10 ' 



B = (z)=A 



(19) 



The coefficient A being a positive quantity, we remark that 
the function l|18|l is normalizable if and only \i B = A> Q. 
This provides a simple derivation of the stability criterion 
that we have used : when A < Q, the fixed point x = is 
stable (the distribution H18|) being not normalizable, the 
stationary distribution is the Dirac function at 0); when 
yl > 0, the fixed point is unstable and the stationary dis- 
tribution is extended. 

The distribution (|18|) is obviously a power law for small 
values of the variable r which is an indication for inter- 
mittency. In figure 13 numerical results (open circles) for 
the probability distribution of the variable r are compared 
with the analytical formula H18|l (dashed line); the coeffi- 
cients A^ B and C have been calculated using Hll|l . The 
agreement between the two curves is excellent as far as the 
power-law behaviour for small values of r is concerned. For 
higher values of r, a discrepancy appears : the assumption 
that the stationary distribution is separable is no more 
valid and a specific analysis for large r is needed |13| . In 
figure^ we plot the probability density of the energy for 
a greater value of the noise : again the small r power law 
is very well described by equation 1)18(1 . 
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Fig. 4. Probability density (in log-log scale) of the variable 



r = vx^ + i^ for 7 = 1, zi = 6 and fc = 1. 
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Fig. 5. Critical amplitude Zi* of the noise versus 7 for a time 
correlation r = 1. The dotted line is a linear fit. 



5 Effect of a colored noise 



In order to study the effect of the temporal correlations 
of the noise, we have performed numerical simulations of 
equation Q), taking the noise to be an Ornstein-Uhlcnbeck 
process. In the stationary regime, the correlation of the 
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Ornstein-Uhlenbeck process is exponential : 

A 



mm)^^exp . (20) 



6 Conclusion 



Using the time series of the simulation results, we estimate ^^-^ ^^^^^^ ^^^^ ^^^^-^^ ^1^^ noise-induced bifurca- 
the Lyapunov exponent (by computing (x/x)) and identify ^-^^ ^ nonlinear oscillator with a fluctuating damping 
the bifurcation threshold. In figure we plot the critical ^^^^ ^^^^ ^^^^^^ ^^^^ instability threshold of the 
value of the noise intensity A for the onset of instability as nonlinear oscillator can not be determined from a stabil- 
a function of 7, the correlation time r being equal to 1 (in -^^ analysis of the moments of the hnearized problem; the 
dimensionless units). We observe that the critical curve ^^^^^^ j^^^^^ ^^^^ ^^^^1^ ^1^^^^ analysis of moments 
has the shape of a straight line. Performing this analysis indicates. This feature, that has already been noticed in 
for different values of r, we have found that the critical ^ ^^^-^^^ ^^^^1^ fmn\ . can be explained by the ex- 
amplitude A, of the noise is almost linear with 7; hence, -^^^^^^ ^f long tails in the probabihty distribution of the 
to a good approximation, the critical curve is given by linearized system. These long tails have a greater influ- 
Z\,(7) = C(t)7. In figure E we plot C{t) as a function ^^^^ ^-^^^^ moments of x which are therefore less and 
of r; again the evolution is well described by a straight jg^g g^^^le. However, these long tails are suppressed for 
hue in the range < r < 1. These resuhs indicate that nonlinear system, which has, therefore, a weU defined 

the structure of the bifurcation diagram for the Ornstein- stability threshold. For the problem considered here, an 
Uhlenbeck noise remains qualitatively the same as that for ^^^^^ calculation of the Lyapunov exponent has allowed 
white noise. ng determine the exact stability diagram of the system. 

Besides, we have shown that a good quantitative descrip- 

lOi ' ' ' ' 1 tion of the system in the vicinity of the bifurcation can 

9- - be obtained if the probability distribution is assumed to 

8- C(x) = 4 + 5.61 X be separable in energy and angle variables. This Ansatz 

allows us to calculate the power law exponent of the sta- 
tionary probability and to predict intermittency. Our the- 
oretical findings agree with numerical simulations. Finally, 
we have studied numerically the effect of temporal corre- 

0.8 1 

lation of the noise by coupling our dynamical system to 
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Fig. 6. Evolution of the critical slope as a function of the an Ornstein-Uhlenbeck process, and have found that the 
correlation time. The dotted line is a linear fit. bifurcation scenario remains qualitatively the same. For 
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this latter problem, exact analj^tical calculations seem to 16. C. W. Gardiner, Handbook of stochastic methods (Springer, 

be out of reach; some quantitative information may, how- 1984) 

ever, be obtained either by using one of the various Marko- 17. H. Riskcn, The Fokker-Planck equation (Springer, 1989) 

vian approximations for colored noise, or by considering ^-^^ Gradstem, I.M. Rhizik, Table of Integrals, Series and 

a special type of random process, such as the dichotomic (Academic Press, 1994) 

Ti • • 1 • I, 1 J • ii J. J. 19. J.P. Bouchaud, A. Comtet, A. Georges, P. Ledoussal, Eu- 

Poisson noise which leads, m the present case, to an ex- ' i o i i 

, , , , , rophys. Lett. 3, 653 (1987) 
actly solvable problem. 

20. E. Ott, J.C. Sommerer, Phys. Lett. A 188, 39 (1994) 
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